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We discuss the detection of entanglement in interacting quantum spin systems. First, thermo- 
dynamic Hamihonian-based witnesses are computed for a general class of one-dimensional spin- 1/2 
models. Second, we introduce optimal bipartite entanglement observables. We show that a bipartite 
entanglement measure can generally be associated to a set of independent two-body spin observables 
whose expectation values can be used to witness entanglement. The number of necessary observables 
is ruled by the symmetries of the model. Illustrative examples are presented. 
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Entanglement is a striking feature of quantum me- 
chanics, revealing the existence of non-local correlations 
among different parts of a quantum system. Entangle- 
ment has been recognized as an essential resource for 
quantum information processing 0. This has provided 
strong motivation for studies probing for the presence 
of naturally available entanglement in interacting spin 
systems 0, 1^- Moreover, the realisation that entan- 
glement can also affect macroscopic properties (such as 
the magnetic susceptibility) of bulk solid-state systems 
0, 0| , has increased the interest in characterizations of 
entanglement in terms of macroscopic thermodynamical 
observables. An observable which can distinguish be- 
tween entangled and separable states in a quantum sys- 
tem is called an entanglement witness |^ . Several differ- 
ent methods for experimental detection of entanglement 
using witness operators have been proposed 0- Entan- 
glement witnesses have recently been obtained in terms of 
expectation values of thermodynamical observables such 
as internal energy and magnetization |^ llCl| , and mag- 
netic susceptibility 01. 

Our aim in this work is two- fold: first, we find 
an entanglement witness for a broad class of interact- 
ing spin-1/2 particles, thus generalizing the result of 
Refs. This is an entanglement witness for all 

spin-1/2 based solid-state quantum computing propos- 
als, such as electron spins in quantum dots ^] and P 
donors in Si 0. While this approach is very general, its 
drawback is that it is sub-optimal, in the sense that it 
does not detect all entangled states. In contrast, in the 
second part of this work, we introduce the concept of an 
optimal bipartite entanglement observable. This allows 
us to construct optimal bipartite-entanglement witnesses 
for qubit systems. The essential idea here is to directly 
relate bipartite entanglement measures and the expecta- 
tion value of spin observables [l3| . 

Hamiltonian-based entanglement witnesses. — An im- 
portant class of spin-based solid-state quantum com- 
puting proposals is approximately governed by diago- 
nal exchange interactions (involving only erf cr" terms, 
where a € {x,y, z} and erf is the Pauli matrix for spin 



i) However, spin-orbit coupling introduces off- 

diagonal terms into the exchange Hamiltonian In 
this case previous results concerning Hamiltonian-based 
entanglement witnesses H, 0, do not apply, since they 
are restricted to the diagonal case. Here we construct an 
appropriately generalized entanglement witness. 

The most general Hamiltonian describing N nearest- 
neighbor coupled spin-1/2 particles in ID is of the 

form H ^ E.Ea,/3G{^,y,^}ffM+i'^>f+i' ^^^^^^ dtf+i = 
(sfi+i)*- There are thus nine independent parame- 
ters for each pair of spins i,i + 1. It is convenient 
to re-express H in terms of a scalar part and sym- 
metric and anti-symmetric parts. In addition we allow 
for the presence of a global external magnetic field B: 

+Yld=i A-(crjXcrj+i)-|-X;,^i(C-crj)(C-(Ti+i), where (T; = 
(erf, fjf , erf ), J a are exchange coupling constants, and we 
assume periodic boundary conditions {ctn+i = ci)- The 
anisotropic term involving A (the Dzyaloshinskii-Moriya 
vector in solid-state physics) typically arises due to spin- 
orbit coupling; |A|/J has been estimated to be in the 
range 0.01 — 0.8 in coupled quantum dots in GaAs [l3| . 
The vector C can arise also due to dipole-dipole coupling 
and other sources. 

We now derive a thermodynamical entanglement wit- 
ness for a system governed by H. Let J — maxcd Jq,|}, 
A = maxa{|v4c|}. Let {X) = Tr(pAr), with p the system 
density matrix. Let u = (H) /N be the per-spin internal 
energy and m = (nix, my, m^) the magnetization vector 



with components = X^i^i i'^t) /-^ ■ Then H yields: 



N\u + B- ml < 



N 
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N 
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^((C.^,)(C-^,+i)> 



Consider an arbitrary separable density matrix p — 
J^kPkPk ^ Pk ® ■■■ ^ Pk ^ where J2kPk = 1 and all 



2 



Pk > 0. It has been shown that for such p, us- 



ing the easily verified facts {cr°'a^_^^) = (af ) ((T^_^_^) 
and J2a=x y z i'^i)'^ — 1' ^^'^ Cauchy-Schwarz 
(CS) inequahty lE.a.^d < (E. a? ^')'/', that 

0- We therefore obtain bounds for the remaining 
two terms. Let Xi = (erf), {xy)i = {<jfa^_^_i), etc. Using 
again the above facts and the CS inequahty we have 



N 



N 



1^ (A • (cT, X (T,+i)} I < ^ I (A • (cr, X I 

N 

<2^^|(yz), + + (a;y),| 

N 

= 2A'^ lUiZi+i + ZiXi+i + XiUi+i I 

i=l 
AT 

< 2Aj2iJ2{aTf)^/^iJ2{a:^,)y/' < 2NA. 



Note that if we assume no symmetry breaking, i.e., 
(af) = (a"), then in fact | ^ti (A • (t. x cr,+i)) | = 0. 

We now obtain an upper bound for the third term. In 
the standard Bloch-sphere parametrization for the indi- 
vidual spin density matrices we have: p = ^(/ -I- n ■ <t), 
where |n| < 1. Then (C ■ cr) = C ■ n < |C| |n| < |C|. 
Therefore, 



N 



N 



\J2iiC- a,){C ■ \<J2\{{C. <T.)(C . <T.+i))| 

i=l i=l 
N 

= Y,\{C-CT,)\\{C-CT,+,)\<N\C\\ 

2=1 

These upper bounds combine to yield the entanglement 
witness: 



W 



B-m|/(J + 2A- 



|C|^). 



(1) 



The numerator consists of macroscopic, observable quan- 
tities. The denominator consists of material parameters. 
We have seen that separability implies W < 1. There- 
fore if l¥ > 1 the system is entangled. When Ja = J 
and the anistropic terms in H are entirely due to the 
spin-orbit interaction, it is possible to relate H to the 
isotropic Heisenberg Hamiltonian via a unitary transfor- 
mation 0, 01 . Applying this transformation to the ex- 
amples of entangled states that are detected by W in 
the case A = C = 0, found in Refs. [1,13) yields exam- 
ples of non-trivial entangled states detected by W when 
A, C 7^ 0. The importance of the witness W is that it is 
directly applicable to a wide class of spin- 1/2 based solid- 
state quantum computing proposals [ill ll^ , where the 



effect of spin-orbit coupling is known to be non-negligible 

M 

Spin-based entanglement witnesses. — Let us turn now 
to the construction of optimal bipartite-entanglement 
witnesses, based on spin observables. Consider a general 
two-body observable R = Y^Rap^siij) (lli i^lj, 

where is a basis for the Hilbert space, i,j enumer- 

ate d -level systems, and a, (3,^, 6 S {0, 1, ...,d— 1}. The 
expectation value of R can generally be written as 



;i?)=5]Tr(RfeV^), 



(2) 



where R(«j) are d x d matrices with elements Rap-fs{ij) 
and p'^ is the two-body reduced density matrix. Eq. Q 
holds for any mixed state and for any d. Here we are 
especially interested in d = 2, i.e., the qubit case. We 
then use the standard basis {|00) , |01) , |10) , |11)} for any 
pair {i,j) of spins, and denote — (0^0^ | p'^ jOiOj), 
Pi2 = (OiOjl/5*-' |Oilj), etc. For an operator R display- 
ing a constant interaction between nearest-neighbor par- 
ticles, the non-vanishing matrices R(zj) are given by 
R(i,i -|- 1) = R (Vi). Moreover, if translation invari- 
ance is assumed in the system then = p. Hence we 

can rewrite Eq. ||2J) as 



;7^) = Tr(Rp), 



(3) 



where TZ = R/N, with N the number of nearest-neighbor 
pairs in the system. 

A convenient bipartite entanglement measure is the 
negativity |l7j . ranging from (no entanglement) to 1 
(maximal entanglement), defined as follows: 



J\f{p) — 2 max(0, — min(/ia)). 



(4) 



where pa are the eigenvalues of the partial transpose 
of the two-particle reduced density operator p, given by 
{aP\ p^^ \^S) = p \ad). We denote the lowest eigen- 
value of p-^-* by p„i which, for composite systems of di- 
mensions 2x2 and 2 x 3jis non-negative if and only if the 
state is separable [l^ll^- Thus, from Eq. |@J, one can 
see that separability implies vanishing negativity 0|. 

The eigenvalue pm, which is the key object of our 
framework, is generally a non-linear function of the ma- 
trix elements of the density operator. However, let us 
first consider the linear case, which occurs for several 
interesting quantum spin systems, as will be illustrated 
later. In this case, we can directly relate p™ to a single 
observable which plays the role of an entanglement wit- 
ness. Indeed, assume that p„i = b=i fa^bPab, where 
fab are constants. Then, defining the matrix elements of 
TZ as Rab — fba we obtain 



p„ = Tr(Rp) = (7e; 



(5) 



Therefore, the observable TZ directly detects the existence 
of bipartite entanglement in the system. (TZ) < im- 



3 



plies bi-partite entanglement, while otherwise the state 
is separable. 

Eq. © can also be established, formally, for those 
cases where depends non-lincarly on p. Indeed, from 
Eq. ||2Jl it follows that 



[U) 



Tr(R'^^p^-*) 



(6) 



where the matrix R"^-^ is defined through (a/3| R"^^ \^5) — 
(7/3|R|ai5), is a unitary matrix which diagonalizes 
p^^ (note that p^-* is Hermitian 18J), and the {pa] de- 
note the four eigenvalues of p^^ . If we choose TZ such 
that 



lyR^^IF^ = diag(ei,e2,e3,e4), 



(7) 



where Cq = 1 for the value of a such that pa ~ A*m and 
Cq, = for the other ones, then as desired (TZ) — pm- 
Hence the expectation value of TZ can be used in general 
as a criterion of separability. However, in the non-linear 
case, if we define TZ through Eq. iQ, TZ itself will be a 
function of the density matrix (since W is) . One would 
then need to measure a complete set of observables to find 
TZ, which just corresponds to quantum state tomography 
0. We show below that in the presence of symmetries 
the number of measurements required to construct TZ can 
be drastically reduced. 

XYZ spin chain in a magnetic field. — In order to 
provide an example of spin-based witnesses in quan- 
tum spin chains, let us consider a parametric family 
of spin Hamiltonians H = ^^-^ i?i.i_|_i(0i, ©i+i), where 

U{Oi) ~ cxp{i6i ■ (Ti), and 



H. 



J u 



(8) 



with periodic boundary conditions assumed, i.e. cr^+i = 
af. Observe that a large class of one-dimensional spin 
models is covered by the Hamiltonian JHl . This family of 
Hamiltonians obeys the constraint \jl,af(jfj^^ = and 
belongs to the subalgebra sm(2)©sm(2) C sm(4). Defining 

— 2 "t+l ^ "i+lj' ^1—2 \"i "i+l ^ "i "i+l)^ 

Zf = ^{a^±a^_^_^), one of the su{2) terms is generated by 
the set of operators {T^ , R~ , Z^} (respectively with co- 
efficients Jx+Jy, Jxy — Jyx, and h, in H) and preserves the 
two dimensional subspace spanned by {|01) , |10)}. The 
other su{2) is generated by {T~ , Rf , Z^} (respectively 
with coefficients Jx — Jy, Jxy + Jyx, and h, in H) and pre- 
serves the other two dimensional subspace spanned by 
{|00),|11)}. The Hamihonian 7Ji,i+i(0i, has the 

same entanglement properties as since they are 

connected through local unitary transformations. There- 
fore, assuming that the initial state has the same sym- 
metry as the Hamiltonian (no symmetry breaking), the 



nearest-neighbor reduced density matrix for an arbitrary 
mixed state reads, in the standard {|00) , |01) , |10) , |11)} 
basis 



a y 

b z 

z* c 

y* d 



(9) 



Positivity of p implies that ad > \y\'^ and be > \z\'^, with 
a,b,c,d > 0. Computing the eigenvalues of p'^^ leads 
to two independent possibilities for the lowest eigenvalue 



,(1) - 



(2) 



(a + d-y^(a-d)2+4|2|')/2, (10) 
'+4\yf)/2. (11) 



4^ 



The condition for entanglement pm < 0, together with 
positivity of p, yields the restrictions 



p^,j' < ad < 6c, 
be < ad. 



p(^) < 



(12) 
(13) 



Note that, since pm is nonlinear in the density matrix 
elements, the resulting witness TZ, obtained from Eqs. © 
and |(7l, will be density matrix-dependent. Indeed, for 



pm , 



TZ = Im(/)i?- + Re(/)T^ 



+ - (I iS> I + (7"" a"") , 
4 



1 (a-d) 



^(a-d)2 +4|z|2 



(14) 



where / = -z/y^{a - dY + 4|z|2, and R^,T+,Z+ were 
defined above. It is seen from this result that the state- 
dependence can be removed if the following constraints 
are obeyed: a — d and z is either real or imaginary. These 
constraints are obeyed, e.g., for the isotropic Heisenberg 
model (see below) . In this case we need to measure just 
one observable in order to determine the entanglement 
properties of the system. Analyzing the second possibil- 
ity, i.e. pr, 



(2) 1,4- • 

Pm , we obtam 



TZ = -Ini(5)i?+ -t- Re(5)T" 



+ -{I ®I - a'' ®(t'') , 



1 



(6-c) 



(15) 



where g = —y/ ^(6 — c)2 + 4|?;|2, and R^,T ,Z were 
defined above. Similarly, TZ is state-independent in 
Eq. H15() for b — e and y either real or imaginary. An 
example of this case is given by the transverse field Ising 
model (see below). But even in the case of the rather 
general Hamiltonian ©, it is clear that instead of full- 
scale quantum state tomography, it suffices to measure 



4 



the elements {a,d,z} or {b,c,y}, in order to construct 
the witness operator TZ. 

Heisenberg model. — Let us consider the constraints 

Jx — Jy — 'Jz — J ^ and Jxy — Jyx — 

/i = in Eq. 0. Then we have the antiferromag- 
netic Heisenberg chain, whose Hamiltonian reads H — 
JT,f=i («+i + f^i^^i+i + ^l^i+i)- H is invariant un- 
der cyclic spin translations and has SU(2) symmetry, 
with H commuting with the total spin components 
^^crf , a e {x,y,z}. The elements of p in Eq. then 
obey further constraints, namely, y = 0, z ~ z* < 0, 
a ~ d = 1 + {crfcrj.^), and b = c = 1 ~ (a^af^^), with 
«crf+i) < li 01(a). It then follows from Eq. ^ 

(2) 

that the eigenvalue /i™ is always non-negative, whence 
entanglement is determined by the eigenvalue iiln , given 
by Eq. (|10|l . Thus, the witness comes from the ob- 
servable in Eq. itTH) . which becomes TZ = {afaf^^ + 
'^i<^i+i + cTiC^i+i + ^)/4 (for any site i). This yields the 
entanglement witness {TZ) = j {u / J + 1) , where we have 
used that, due to the translation invariance and SU{2) 
symmetry, the correlation functions satisfy the relations: 
(afaf+i) = (afaf+i) = (afcTf+i) = u/i3J) :i3](a). 
Remarkably, this spin-based witness is (upto an irrele- 
vant prefactor) precisely the Hamiltonian-based entan- 
glement witness found in Ref. 8] (and in our general- 
ized Hamiltonian-based result above). Since our spin- 
based approach is optimal for bi-partite entanglement 
{{TZ) is essentially the negativity), this is a proof that 
the Hamiltonian-based witness Q detects all bi-partite 
entangled states. 

Transverse Field Ising model. — As a final example we 
analyze the ferromagnetic one-dimensional Ising chain 
in the presence of a transverse magnetic field. This 
model corresponds to taking = —XJ and h — — J in 
Eq. ((HJ, with J > and all the other couplings vanishing: 
H = -JEti Then y = y* {Z, sym- 

metry of H) and b = c (translation symmetry of H) [16| . 
From the analysis of the thermal correlation functions, 
which can be obtained analytically j2oj . it can be shown 
that, in contrast to the Heisenberg case, entanglement 
is now determined by the eigenvalue /im^ in Eq. 
From Eq. 115|l . our spin-based entanglement observable 
is then TZ = - jicrfaf^^ - crfcrf+i + erf crf+i - 1) for any 
site i. The expectation value of this observable can be 
determined from measurements of the nearest-neighbor 
spin-spin correlations {afaf^i). This can be done, e.g., 
by inelastic neutron scattering We note that it was 
shown in Ref. IBj that spin-spin correlation functions can 
act as entanglement witnesses in bulk solids. Our witness 
operator TZ is optimal, so that it can, moreover, yield 
precise macroscopic predictions. For instance, we plot in 
Fig. n the witness as a function of /3 = 1/kT for sev- 
eral values of A, where k is the Boltzmann constant and 
T is the temperature. From this figure, we can obtain 
the exact critical temperature Tc for the entanglement- 



separability transition. For A = 1, we have pc ~ 1-93 or 
kTc « 0.51 (in units such that J = 1). This temperature 
can be compared to the value 0.41 obtained in Ref. 
by using an energy witness. The reason for the small 
difference is that the energy witness of Ref. is derived 
from an entanglement bound and, despite being a good 
approximation, neglects some entangled states which are 
detected by TZ. 




FIG. 1: Witness for the transverse field Ising model. The 
entanglement-separability transition temperature is indicated 
by the intersection with the horizontal axis. 
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